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Abstract 



In this paper we present a unifying geometric framework for mod- 
eling various sorts of physical network dynamics as port-Hamiltonian 
systems. Basic idea is to associate with the incidence matrix of the 
graph a Dirac structure relating the flow and effort variables asso- 
ciated to the edges, internal vertices, and boundary vertices of the 
K*" I graph. This Dirac structure captures the basic conservation/balance 

^y ■ laws of the system. Examples from different origins such as consensus 

^—v I algorithms and coordination control strategies for multi-agent systems 

p^ ■ share the same structure. The framework is extended to fc-complexes 

primarily motivated by the discretization of continuous conservation 
laws. 



1 Introduction 



/\ ' Discrete topological structures arise abundantly in physical systems mod- 

j^ ■ eling. A classical approach to modeling and analysis of electrical circuits, 

dating back to Kirchhoff, is through the circuit graph. Similar approaches 
apply to many other cases, including e.g. mass-spring-damper mechanical 
systems, multi-body systems, hydraulic networks, chemical reaction net- 
works, and power systems. Common feature of all these cases is that the 
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discrete structures, in particular graphs, are blended with dynamical rela- 
tions, leading to various sorts of network dynamics. 

During the last two decades network dynamics has received ever-increasing 
attention, with inputs, among others, from graph theory, multi-agent sys- 
tems, dynamical systems, as well as statistical mechanics. In this paper we 
wish to lay down a general geometric framework for defining physical dynam- 
ics on directed graphs. We will restrict attention to fixed directed graphs, 
that is, we do not study the possible dynamics of the graphs themselves, as 
e.g. studied in random graph theory or statistical mechanics. Furthermore, 
our graphs physically correspond to some sort of conservation (or balance 
laws), relating real- valued variables associated to the vertices to those cor- 
responding to the edges of the graph. Dynamics arises because of storage 
(accumulation) of quantities corresponding to some of these variables, ei- 
ther at the vertices or at the edges (or both). Finally, in order to capture 
interaction with the environment and interconnection to other networks we 
explicitly consider boundary vertices, corresponding to the terminals of the 
network. 

The generalized Hamiltonian nature of the resulting dynamical models 
is due to the common assumption that the flux variables corresponding to 
storage at the vertices and/or edges are derivable from some energy function 
(defining the energy corresponding to storage) . Next to these energy- storage 
constitutive relations, the remaining internal variables will be related by 
static energy- dissipating constitutive relations. This will imply that the total 
stored energy (the Hamiltonian function) itself satisfies a conservation law: 
the increase of the total energy is equal to externally supplied power (through 
the boundary vertices of the graph), minus the power lost in the dissipative 
elements (corresponding to some of the edges or vertices of the graph). The 
resulting systems fall within the class of port-Hamiltonian systems., as coined 
and explored in e.g. [3511321 [36l ll3j. They are truly ^enerafeec? Hamiltonian 
systems, since they allow for energy-dissipation and interaction with the 
environment. 

From a geometric point of view the generalized Hamiltonian structure 
of the network dynamics is defined, apart from its Hamiltonian function 
and energy-dissipative relations, by a Dirac structure. This Dirac structure 
(generalizing the symplectic or Poisson structure from classical mechanics) 
is directly defined by the incidence matrix of the directed graph, and thus 
captures the conservation laws. In fact, we will show how a directed graph 
gives rise to three canonically defined Dirac structures on its vertex and edge 
spaces. The first two of them only differ in the different role of the boundary 
vertices, while the third, the Kirchhoff-Dirac structure, captures the special 



case where no storage or dissipation is associated with the vertices of the 
graph (corresponding to Kirchhoff's current laws). 

We win ihustrate this framework on some of the physical examples men- 
tioned above. Furthermore, we will show how the same port-Hamiltonian 
structure is shared by network dynamics with a different origin, such as 
consensus and coordination control algorithms. 

While perhaps all examples given in the paper are relatively simple, and 
can be approached from other angles as well, we believe that a major contri- 
bution of the paper resides in the identification of a common mathematical 
structure in all these examples, which is moreover closely related to classical 
Hamiltonian and gradient dynamical systems. 

In the last part of the paper we will sketch how the geometric framework 
as developed for directed graphs can be extended to arbitrary /c-complexes. 
By considering standard simplicial complexes this allows a structure-preserving 
spatial discretization of distributed-parameter physical systems, otherwise 
described by partial differential equations. This is illustrated on the 2-D 
Maxwell equations and a general form of discretized diffusive systems. Fur- 
thermore, we explore its potential fro constructing generalized consensus 
algorithms. 

Preliminary work regarding Sections 3.4 and 3.5 can be found in [39], 
and regarding Sections 5.1, 5.2 and 6.2 in [371 [38] • 



2 From directed graphs to Dirac structures 

In this section we will define two, closely related, canonical Dirac structures 
on the combination of the vertex, edge and boundary spaces of a directed 
graph, and their dual spaces. Before doing so we first recall some basic 
notions of graph theory and Dirac structures. 

2.1 Directed graphs and their vertex and edge spaces 

First we recall some standard definitions regarding directed graphs, as can 
be found e.g. in [5j. A directed graphic Q consists of a finite set V of vertices 
and a finite set £ of directed edges, together with a mapping from 8 to 
the set of ordered pairs of V, where no self-loops are allowed. Thus to any 
branch e ^ £ there corresponds an ordered pair (f;, u;) € V x V (with v ^ w), 
representing the tail vertex v and the head vertex w of this edge. 



'^ Sometimes called a multi- graph since we allow for the existence of multiple branches 
between the same pair of vertices. 



A directed graph is completely specified by its incidence matrix B, which 
is an A^ X M matrix, N being the number of vertices and M being the number 
of edges, with (i, j)-th element equal to —1 if the j-th edge is an edge towards 
vertex i, equal to 1 if the j-th. edge is an edge originating from vertex i, and 
otherwise. Since we will only consider directed graphs in this paper 'graph' 
will throughout mean 'directed graph' in the sequel. 

Given a graph, we define its vertex space Aq as the vector space of all 
functions from V to some linear space TZ. In the examples, TZ will be mostly 
TZ = R OT TZ = W^. In the first case, Aq can be identified with M . Further- 
more, we define its edge space Ai as the vector space of all functions from 
£ to the samqj linear space TZ. Again, if 7^ = M then Ai can be identified 
with M^'-^. 

The dual spaces of Aq and Ai will be denoted by A'^, respectively A^. 

Remark 2.1. Hence A^ can be represented as the linear space of functions 
from V to the dual linear space TZ* , and h^ as the linear space of functions 
from £ to TZ* . Thus the duality pairing between f £ Aq and e £ A^ is given 
as 

< f\e>=Y,< fiv) I e{v) > , 
vev 

where < \> on the right-hand side denotes the duality pairing between TZ 
and TZ* , and a similar expression holds for / G Ai and e € A^. 

The incidence matrix B of the graph induces a linear map B from the 
edge space to the vertex space as follows. Define i? : Ai — >■ Aq as the linear 
map with matrix representation 

B^I 

where I : TZ ^ TZis the identity map and CS) denotes the Kronecker product. 
B will be called the incidence operator. For 7^ = R^ the incidence operator 
is given in matrix form as B CS) I3, with I3 the 3 x 3-identity matrix, while 
for 7^ = M the incidence operator reduces to the linear map given by the 
matrix 13 itself. In this case {TZ = M) we will throughout use B both for the 
incidence matrix and for the incidence operator. 
The adjoint map of B is denoted as 

B* : A°^ A\ 



^In principle we could also associate with the edges a linear space TZ' which is different 
from the space 7?. associated with the vertices. In that case the definition of the incidence 
operator needs an additional linear map from 7?,' to 7?.. 



and is called the co-incidence operator. For 7^ = M'^ the co-incidence opera- 
tor is given by B^ (8> I3, while for 7^ = M the co-incidence operator is simply 
given by the transposed matrix B . In the latter case we will throughout 
use B^ both for the co-incidence matrix and the for co-incidence operator. 
We will use the terminologjo flows for the elements of Aq and Ai (nota- 
tion /o and /i), and efforts for the elements of their dual spaces A*^ and A^ 
(notation e^ , respectively e^). 

2.2 Open graphs 

An open graph Q is obtained from an ordinary graph with set of vertices V 
by identifying a subset V{, C V of N^^ boundary vertices. The interpretation 
of Vb is that these are the vertices that are open to interconnection (i.e., 
with other open graphs). The remaining subset Vi := V — Vb are the Ni 
internal vertices of the open graph. 

The splitting of the vertices into internal and boundary vertices induces 
a splitting of the vertex space and its dual, given as 

Ao = AoieAofe 
A" = AO^eAO'' 

where Aoi is the vertex space corresponding to the internal vertices and Agfe 
the vertex space corresponding to the boundary vertices. Consequently, the 
incidence operator B : Ki ^ Kq splits as 

B = Bi®Bb 

with i?j : Ai — )■ Aoj and Bb : hi ^ Aq;,- 

Furthermore, we will define the boundary space Af, as the linear space of 
all functions from the set of boundary vertices Vb to the sameQ linear space 
TZ. Note that the boundary space A^ is isomorphic to the linear space Agb, 
and that using this isomorphism the linear mapping Bb can be also regarded 
as a mapping 

Bb-.Ki^ Kb 



^This terminology steins from port-based and bond-graph modeling |26| . where it has 
a shghtly more specific connotation then in our case. The space Ao is also called the space 
of 0-chains, while the elements of Ai are called the 1-chains. Furthermore, the dual spaces 
A° and A^ are called the space of 0-cochains, respectively 1-cochains. 

*Again, in principle we could also associate with the boundary vertices a linear space 
which is different from the space TZ associated with the vertices and the edges. 



called the boundary incidence operator. Nevertheless, we will be careful in 
distinguishing the two isomorphic linear spaces A^ and Ag;, because of their 
different interpretations in physical examples. The dual space of A^ will be 
denoted as A^. The elements /^ € A^ are called the boundary flows and the 
elements e^ G A^ the boundary efforts. 

Remark 2.2. Another interesting case, which for clarity of exposition will 
not be considered in this paper, is to allow for interaction between graphs 
through boundary edges. 

2.3 Dirac structures 

Recall ( [351 |8l |32] ) the definition of a (constaniO) Dirac structure. Consider 
a vector space T with dual space T*. As before, the variables f € T are 
called the flow variables, while the conjugate variables e E J-"* are called 
the effort variables. Define on the total space J-" x J^* the indefinite inner 
product ^ •, • 3> as 

< {fa,ea),{fb,eb) »:=< ea \ ft > + < et \ fa >,fa,fb e J^,ea,eb G T* 
where < • | • > denotes the duality product between J^ and J-"*. 
Definition 2.3. A subspace T> <Z T x T* is a Dirac structure if 

V = V^ (1) 

where -^ denotes the orthogonal complement with respect to <C •,• ^. 

In the finite-dimensional case an equivalent, and often easier, character- 
ization of Dirac structures is given as follows (see e.g. [3 [13] for a proof). 

Proposition 2.4. A subspace T> d T y. T* is a Dirac structure if and only 
if the following two conditions are satisfied: 

«<e|/>=0, forall{f,e)eV 
(ii) dim D = dim J^ 



^This definition can be extended [12118] to (non-constant) Dirac structures on manifolds: 
a Dirac structure I? on a manifold M is defined as a vector subbundle of the Whitney 
sum TM ffi T*M such that for each m £ M the linear space D{m) C TmM x T^M 
is a constant Dirac structure. This will be needed, for example, in the treatment of 
mechanisms in Section 3.3. 



Note that the first equation in ([2|) can be regarded as a power- conservation 
property. Indeed, if / and e are power-conjugated variables, that is, < e | / > 
denotes power, then the first equation of ([2]) states that the total power in 
the system is equal to zero. 

A special type of Dirac structure is defined as follows: 

Definition 2.5. A Dirac structure T) <Z F x T* is separable if 

<ea\ fb>=0, for all {fa, e^), {fb, Cb) G V (3) 

Note that the definition of a Dirac structure already implies the weaker 
property 

< Ca I /fe > + < eb I /a >= , for all (/„, €«), {fb, Cf,) G V (4) 

Separable Dirac structures have the following simple geometric characteri- 
zation. 

Proposition 2.6. Consider a separable Dirac structure T> d J- x T* . Then 

V = JCxJC^ (5) 

for some subspace K <Z T, where K^ = {e G J^* |< e | / >= 0,V/ G /C}. 
Conversely, any subspace T> as in ^ for some subspace K d J- is a separable 
Dirac structure. 

Proof. It is immediately seen that any subspace K, x K,^ satisfies ([3]), 
and is a Dirac structure since it satisfies a fortiori @, and has dimension 
equal to dim F. Conversely, let the Dirac structure V satisfy ^ . Define the 
following subspaces 

J-Q = {/ G ^ |(/, 0) G P} -Fi = {/ G -F I 3e G F* s.t. (/, e) G V} 

^0 = {e G .F* I (0, e) G P} ^i = {e G F* [ 3/ G F s.t. (/, e) G V} 

It is readily seen [9] that for any Dirac structure 8i = {Tq)^,£q = (Fi)^. 
We will now show that ([3]) implies that Fq = Fi =: /C (and hence £q = £i =: 
K,^). Clearly, Fq d Ti. Let now {fa, Go) G V and thus /a G F"!. Then for 
ah {fb,eb)^V 

< (/a,0), {fh,eh) >:=< ea I /fe > + < I /a >=< Ca I /fe >= 

by (l3|). Hence, also {fa, 0) G P and thus fa G Fq. By definition Fq x £"9 C T), 
and hence K. x K,^ C P. Finally, since the dimension of /C x K,^ equals the 
dimension of F equality results. ■ 



Note that ^ can be regarded as a generalized statement of Tellegen's 
theorem for electrical circuits (with / denoting the vector of currents, and 
e denoting the vector of voltages). 

A typical form of a separable Dirac structure, which will be frequently 
used in the remainder, is the following. 

Proposition 2.7. Let ^ : V — )• W 6e a linear map between the linear spaces 
V and W with adjoint mapping A* : W* — )• V* , that is 

<w*\Av >=< A*w* \v> (6) 

for all V ^ V,w* G W* (where, as before, < • | • > denotes the duality 
product between the dual spaces W and W* , respectively V and V*). Identify 
(V X W)* = V* xW*. Then 

V := {{v,w,v*,w*)^{VxW)x{V* xW*)\ 
Av = w,v* = -A*w*} 

is a separable Dirac structure. 

Proof. Define /C:= {{v,w) eVxW\Av = w}. Then /C^ = {{v*,w*) € 
V* xW* \v* = -A*w*}. m 

A key feature of Dirac structures is that their composition is again a 
Dirac structure (in contrast with symplectic or Poisson structures, where 
this is not generally the case). Let Va C Ta x FcX F\ x T* and Vb C 
^B X FcX F^ X F* be two Dirac structures with partially shared space of 
flow and effort variables Fc, respectively F*. Define their composition as 

VaoVb = {{fA, CA, /b, es) (^FaxFbxF\xF*b\ 3(/, e) € 7"^ x F* s.t. 

if A, CA, /, e) E Va, ifs, es, -/, e) G Va} 

(8) 

It has been shown in [71 [3T] that the composition Da o Db of any two Dirac 
structures is again a Dirac structure. Separable Dirac structures have the 
same compositional property: 

Proposition 2.8. Let Da C FaxFcXF^x F* and Db C FbxFcXF^x F* 
be two separable Dirac structures given as 

Di = JC^xlCi ,i = A, B, 



where ICi C J-i x Tc^i = A,B. Define the composition 

ICaoJCb = {{fAjB) e TaxTb \3feTc s.t. UaJ) e /C^, (/b, -/) e ICb} 
Then the composition Da ° T^B is the separable Dirac structure 

VaoVb = (ICa o ICb) x (/Ca o ICb)^ (9) 

For explicit equational representations of compositions of Dirac struc- 
tures we refer to [7]- 

The compositionality property of Dirac structures is a key ingredient of 
port-Hamiltonian systems theory, which implies that the standard intercon- 
nection of port-Hamiltonian systems results in another port-Hamiltonian 
system with Dirac structure being the composition of the Dirac structures 
of the component port-Hamiltonian systems, and Hamiltonian equal to the 
sum of the Hamiltonians of the component systems. 

2.4 The graph Dirac structures 

We now have all ingredients to define Dirac structures corresponding to the 
incidence structure of a directed graph. 

Definition 2.9. Consider an open graph Q with vertex, edge and boundary 
spaces, incidence operator B and boundary incidence operator Bh ■ The flow- 
continuous^ graph Dirac structure I)f{Q) is the following subspace of all flow 
and effort variables 

^fiS) ■■= {ifl,e\fo^,e^\fb,e'')e 

Ai X Ai X Aoj X A0» X Ab X A^ I (10) 

B^fl = foi, Bbfi = fb, el = -5*6°^ - i?*e^} 

The effort- continuous graph Dirac structure T)(,{Q) is 

-DeiG) := {(/i,ei,/o,eO,/,,e^)G 

Ai X A^ X Ao X A° X Afe X A^ I (11) 

Bifi = /o, Bkfi = fob + fb, e' = -B*e\ e' = 6°^} 



The terminology flow- continuous and effort- continuous derives from the fact that in 
the first case the boundary flows /j, are exclusively linked to the edge flows /i, while in 
the second case the boundary efforts e^ are determined by part of the internal vertex 
efforts e". Note that the space of involved flow and effort variables for VfiQ) and I>e(t/) 
is different. 



It directly follows from Proposition 12.71 that both 'Df{Q) and Vg^Q) are 
separable Dirac structures. Note that T)f{Q) and ^^{Q) only differ in the 
role of the boundary flows and efforts, and that 'Df{Q) = T)(,{Q) if there are 
no boundary vertices. 

2.5 Interconnection of open graphs and composition of graph 
Dirac structures 

Interconnection of two open graphs Q°^ and Q" is done by identifying some 
of their boundary vertices, and equating (up to a minus sign) the boundary 
efforts and flows corresponding to these boundary vertices, resulting in a 
new graph. For simplicity of exposition consider the case that the open 
graphs have all their boundary vertices in common, resulting in a (closed) 
graph with set of vertices Vf U Vf U V, where V ■=Vt^ = V^ denotes the set 
of boundary vertices of both graphs. 

The incidence operator of the interconnected (closed) graph is obtained 
as follows. For simplicity of notation consider the case that 7^ = M. Let Q^ 
have incidence operators 



B^ 



Bi 



j = a,(] 



The incidence operator B of the interconnected graph is then given as 



B 



'Bf 


" 





Bf 


B'i 


Bt\ 



(12) 



corresponding to the interconnection constraints on the boundary potentials 
and currents given by 



J>a 



„6/3 



fb+f^ 







(13) 



Of course, several extensions are possible. For example, one may retain the 

Vf^ as being boundary vertices 



set of shared boundary vertices Vfe := 

(instead of internal vertices as above) by extending (flSl) to 



V« 



Jba 



.bP 



e\ f^ + f[: + fb = 0, 



(14) 



with ffj, e* the boundary flows and efforts of the interconnected graph. 

Comparing the interconnection of open graphs with the composition of 
their graph Dirac structures (see e.g. Proposition 12. Sh it is readily seen that 



the flow/effort-continuous graph Dirac structure of an interconnected graph 
equals the composition of the flow/effort-continuous graph Dirac structures 
of ^" and Q^; we leave the straightforward proof to the reader. 

3 Port-Hamiltonian systems on graphs 

In the first subsection we will describe how port-Hamiltonian dynamics can 
be defined with respect to the canonical graph Dirac structures defined 
above. In the subsequent subsections this will be illustrated on a num- 
ber of typical examples, ranging from mass-spring systems and mechanisms 
to consensus algorithms. 

3.1 Definition of port-Hamiltonian systems with regard to 
the graph Dirac structures 

In this subsection we will apply the general definition of port-Hamiltonian 
systems with regard to an arbitrary Dirac structure, see e.g. [Ml [32], to the 
canonical graph Dirac structures as defined above. 

Consider for clarity of exposition throughout the effort-continuous graph 
Dirac structure 'De{Q) involving the flow and effort variables 

ifi,e\fo,e^fb,e') G Ai x A^ x Aq x A" x A^ x A'' 

(The exposition is directly repeated for the flow-continuous graph Dirac 
structure T)f{Q).) A port-Hamiltonian dynamics is specifled by deflning be- 
tween all the internal conjugate flow and effort variables (/i, e^, /o, e'^) either 
an energy- storing relation, or a purely dissipative relation. In general, an 
energy-storing relation between a vector of flow variables / and a conjugate 
vector of effort variables e is of the form 



or dually 









where x is a vector of energy variables (of the same dimension as / and e) , 
and H{x) is any function, representing the energy stored in the system. 

Furthermore, a dissipative relation between a vector of flow variables / 
and a conjugate vector of effort variables e is any static relation 

i?(/,e) = 



satisfying < e | — / >> for all (/, e) satisfying R{f,e) = 0. Often the 
dissipative relation will be of an input-output type, that is, / = —R{e) or 
dually e = -R{f). 

Thus a port-Haniiltonian dynamics is defined by adding to the linear 
relations imposed by the graph Dirac structure constitutive relations be- 
tween all the internal effort and flow variables, either of energy-storing or 
of dissipative typqj It is clear that this leaves many possibilities for defining 
port-Hamiltonian dynamics. In particular, energy-storage, respectively dis- 
sipation, can be associated to the vertices or to the edges, or to both. The 
examples presented in the next subsections cover a number of these different 
possibilities. 

The interpretation of the flow/effort-continuous graph Dirac structure 
as describing discrete conservation or balance laws becomes more clear from 
the above description of port-Hamiltonian dynamics. For example, consider 
for the effort-continuous graph Dirac structure the case of energy storage 
associated to all the edges as well as to the all the vertices: 

x' = e\ h = -^ix') 

xo = -/o, 6° = g (xo) 

for state variables x^ € A^ and xq € Aq, and energy functions H^ and Hq. 
Then the relations imposed by the effort-continuous graph Dirac structure 
imply 

Xo + Bifi = 0, ±^ + B*e^ = 

expressing discrete conservation (or balance) laws between the storage of the 
quantities xq associated to the vertices and the flow /i through the edges, 
respectively between the storage of the quantities x^ associated to the edges 
and the effort e^ at the vertices. The mass-spring system discussed in the 
next subsection will be of this type. 

Finally, we note a fundamental property of any port-Hamiltonian dy- 
namics. Let H{x) denote the total energy of the port-Hamiltonian system. 
Then because of the power-conserving property of the Dirac structure, and 
denoting the flows and efforts of the dissipative elements by /ij, e^, 

j^H{x) =< ^{x) \x>=<e''\fR> + <e'\h ><< e' \ f, > (15) 



^Hence port-Hamiltonian dynamics generalizes both classical Hamiltonian dynamics 
(with no energy-dissipation), as well as gradient systems (where there is in general no 
oscillation between different energies and energy-dissipation does take place); see [34) and 
the references quoted therein. 



Hence the total energy itself satisfies a conservation law: its increase is 
equal to the externally supplied power < e^ \ fb > minus the dissipated 
power — < e \ fji >. Note that individual terms in the summation 
< e^ \ fb >= ^^ e^{vb)fb{vb) (summation over all boundary vertices Vh) 
have dimension of power but do not necessarily correspond to power; see 
the discussion in |32l |3l] . 



3.2 Mass-spring systems 

The first way of modeling a mass-spring system as a port-Hamiltonian sys- 
tem on a graph Q is to associate the masses to the vertices, and the springs 
to the edges of the graph. (Alternatives will be discussed later on.) 

Thus let us consider a graph Q with N vertices (masses) and M edges 
(springs), specified by an incidence operator B. First consider the situation 
that the mass-spring system is located in one-dimensional space 7^ = R, and 
the springs are scalar. A vector in the vertex space Aq then corresponds to 
the vector p of the scalar momenta of all N masses, i.e., p G Aq = M^. 
Furthermore, a vector in the dual edge space A^ will correspond to the total 
vector q of elongations of all M springs, i.e., g G A^ = M . Next ingredient 
is the definition of the Hamiltonian (stored energy) // : A^ x Aq — )• M, 
which splits for a mass-spring system into a sum of the kinetic and potential 
energies of each mass and spring 



where p^ is the momentum corresponding to the mass associated to vertex 
V and Qe is the elongation corresponding to the spatial spring associated to 
edge e. The summations are over all vertices u € V and all edges e £ £. 
Furthermore, m„ are the mass parameters and Vg are the potential energies 
of the springs (e.g., in the case of a linear spring Ve{qe) = \keq1, with k^. the 
spring constant of the linear spring corresponding to edge e) . 

In the absence of boundary vertices the dynamics of the mass-spring 
system is then described as the port-Hamiltonian system 




-B 



B^- 




dR 



di 
dl 
dp 



{q,p) 



defined with respect to the graph Dirac structure T>e{Q) 
that in fact the skew-symmetric matrix 



J:-- 




-B 



(16) 
Vf{g). Note 



(17) 



defines a Poisson structure on the state space A^ x Aq. 

The inclusion of boundary vertices, and thereby of external interaction, 
can be done in different ways. The first option is to associate boundary 
masses to the boundary vertices. Considering the effort-continuous graph 
Dirac structure I'e(^) we are then led to the port-Hamiltonian system 



P = -B^iq,p) + 
e' = [0 l]^iQ,P) 



(18) 



where we have reordered the components of the momentum vector p in 
such a way that its last components correspond to the boundary masses. 
Here /^ € A^ are the external forces exerted (by the environment) on the 
boundary masses, and e^ € A* are the velocities of these boundary masses. 
Another possibility is to start from the fiow-continuous graph Dirac 
structure 'Df{Q). In this case there are no masses associated to the bound- 
ary vertices, and we obtain the port-Hamiltonian system (with pi the vector 
of momenta corresponding to the masses associated to the internal vertices, 
and Hamiltonian H{q,pi)) 



q = Bjmq^p.)+Bl 



SfA'diFtJ^^b 



e 



Pi = -Bi^{q,pi) (19) 

h = B^^iq,p.) 

with e G A the velocities of the massless boundary vertices, and /;, G A^ 
the forces at the boundary vertices as experienced by the environment. 

Note that in this latter case the external velocities e^ of the boundary 
vertices can be considered to be inputs to the system and the forces fb to 
be outputs; in contrast to the previously considered case (boundary vertices 
corresponding to boundary masses), where the forces fb are inputs and the 
velocities e^ the outputs of the systenif|. 

The above formulation of mass-spring systems in 7^ = M directly extends 
to 7^ = M^ by using the incidence operator B = 13 ® I3 and its adjoint as 
defined before. 



*0f course, one could also consider the case where some of the boundary vertices are 
associated to masses while the remaining ones are massless. This will lead to a hybrid 
representation where the inputs consist of a part of the boundary forces and the comple- 
mentary part of the boundary velocities (and the outputs are all the remaining variables). 



Finally, we remark that in the above treatment we have considered 
springs with arbitrary elongation vectors g G A^. For ordinary, straight, 
springs the vector q of elongations is given as, q = B Qc, where qc G A*^ de- 
notes the vector of positions of the vertices. Hence in this case q G imi?^ C 
A^. Note that the subspace iraB^ x Aq C A^ x Aq is an invariant subspace 
with regard to the dynamics (jlSp or ()19p . We will return to this issue in 
Subsection 14.11 

3.2.1 Mass-damper systems 

The above set-up immediately extends to mass-spring-damper systems by 
associating dampers to some of the edges and springs to the remaining ones. 
For brevity we only consider the case of a pure mass- damper system, (no 
springs), with massless boundary vertices. In this case the flow-continuous 
graph Dirac structure yields the following equations 

B^fi = -p 

Bbfi = fb (20) 

e' = -Bl%{p)-Ble^ 

where /i, e^ are the flows and efforts corresponding to the dampers (damping 
forces, respectively, velocities). For example, for linear dampers the consti- 
tutive relations are given as /i = —Re^, where R is a positive semi-definite 
diagonal matrix. Substitution into (j20p then yields the port-Hamiltonian 
systenu 

p = -B,RBl%{p)-B.RBle^ 
h = B,RBj%{p) + BlRBle' 

where, as before, the inputs e^ are the boundary velocities and /{, are the 
forces as experienced at the massless boundary vertices. 
Note that the matrix 



C:-- 



Bi 
Bb 



R[BI Bl] 



^Recall again the generalized Hamiltonian nature of port-Hamiltonian systems. In this 
case the system has only one type of energy storage, i.e., kinetic energy, and thus no 
oscillatory behavior is present. In fact, the mass-damper system can be also represented 
as a gradient dynamical system with respect to an inner product defined by the mass 
parameters and a potential function defined by the damping parameters; see also [34] for 
further information on the relation between port-Hamiltonian and gradient systems. 



is the weighted Laplacian matrix of the graph Q (with weights given by the 
diagonal elements of R) . It is well-known [5] that for a connected graph the 
matrix C has exactly one eigenvalue 0, with eigenvector 1, while all other 
eigenvalues are positive. This implies specific properties of the mass-damper 
dynamics, to which we will return in the context of the consensus algorithm 
(Subsection [37 



3.3 Three-dimensional mechanisms 

This section presents an extension of mass-spring-damper systems in TZ or 
TZ^ to spatial mechanisms, that is, networks of rigid bodies related by joints. 
In this case, the linear space 7?. is given by TZ := se*(3), the dual of the Lie 
algebra of the Lie group SE{3) describing the position of a rigid body in 

A mechanism (or multibody system) is a mechanical system consisting 
of rigid bodies related by joints (defined as kinematic pairs) restricting the 
relative motion between the rigid bodies. The reader may find numerous 
references about their definition and analysis in [2] [29], using different ge- 
ometric representations of rigid body displacements. In this paper however 
we shall follow closely the exposition in |22j|15j . which is based on the Lie 
group of isometrics in Euclidean space M^. 

The basic topology of the mechanism is described by a directed graph, 
called the primary graph, whose vertices correspond to the rigid bodies and 
whose edges are associated with the kinematic pairs. This is similar to the 
mass-spring or mass-damper systems described in Section 13.21 with the dif- 
ference that the dynamical system associated with each vertex corresponds 
to rigid body dynamics instead of point-mass dynamics, and that the edges 
are in first instance associated with kinematic constraints between the bod- 
ies instead of springs or dampers. We shall see how (spatial) springs may 
be included in second instance. 

3.3.1 The rigid body element 

The configuration space of a rigid body is the Lie group of isometrics in 
Euclidean space M^, called the Special Euclidean Group and denoted by 
SE{3) 3 Q (also called the space of rigid body displacements). Using the 
momentum map associated with the action of SE (3) on its cotangent bundle 
T*SE{3), following for instance [El chap. 4], one may define the state 
space of the rigid body as SE (3) x se* (3) 3 {Q, P) by means of the left 
trivialization, where P is called the momentum in body frame. 



The kinetic energy of a rigid body is defined by 

K{P) = ^{P,[l'y\p)) (22) 

wliere / : se (3) — >■ se* (3) is is a symmetric, positive-symmetric isomor- 
phism, called the inertia operator of the rigid body in the body frame.The 
potential energy of the rigid body is defined by a function U{Q) of the dis- 
placement Q. The potential energy may be due to gravity or may be zero 
in the case of the Euler-Poinsot problem. 

We assume that the rigid body is subject to an external force which is 
expressed as an element We G se* (3), called force in fixed frame [16jJ (also 
called wrench in fixed frame [15j) and obtained by the right trivialization 
of T* SE (3). We shall associate a conjugate velocity to this external force, 
the velocity of the body T^ in fixed frame [16] (also called the twist in fixed 
frame [E]), and obtained by the right trivialization of TSE (3). 

The dynamical equations of the rigid body elements may then be written 
as a port-Hamiltonian system [35] |19^ eqn. (1.37)]: 



d 



Te = ( AdQ 




^t[p )-{ -t*Lq -px ) \ {Py\p) )^\Adl ' ^^ 



(23) 



where TLq denotes the tangent map to the left translation (mapping the 
velocities T G se (3) in body frame into the velocities v € TqSE (3)), T*Lq 
denotes its dual map (mapping forces F G TqSE (3) into forces in body 
frame W G se* (3)), Adg denotes the adjoint representation (mapping ve- 
locities in body frame into velocities in fixed frame), Adp denotes the adjoint 
map to Adq , while finally x is defined by the coadjoint representation of the 
Lie algebra se (3) as follows 

W xT = ad*rpW 

for any {W, T) G se* (3)xse(3). The Dirac structure of the port-Hamiltonian 
system (|23|) is thus specified as 



'DRBiQ) = {{v,W,T,,F,T,We)G 

TqSE{3) X se*(3) x se(3) x T^SE{3) x se(3) x se*(3)s.t. 

W ) = [ -T*Lo -pl) [t )^[Ad*Q) ^'=' 



Te = {0 ^dg) (^)} 



(24) 

(Note that this is a non-constant Dirac structure |121[8l[9] on SE{3).) 

In this way we have associated with every vertex of the primary graph 
of the mechanism the pair {We, T^) G se* (3) x se (3) (i.e., the wrench and 
twist of the rigid body in fixed frame), and the dynamical system (I23p . 

3.3.2 The kinematic pair 

Constraints between the rigid bodies of the mechanism will be specified 
by kinematic pairs corresponding to each edge of the primary graph. A 
kinematic pair is the idealization of a set of contacts that occur between 
two rigid bodies at some configuration of the bodies. It constrains the 
possible relative twists between the bodies as well as the possible transmitted 
wrenches. The wrench W transmitted by a kinematic pair is constrained 
to a linear subspace of the space of wrenches se*(3) called the space of 
constraint wrenches and denoted by CW. A relative twist between the two 
bodies is allowed by the kinematic pair when it produces no work with 
any transmissible wrench. The relative twist is thus constrained to a linear 
subspace FT of the space of twists se(3), called the space of freedom twists. 
Since an ideal kinematic pair is workless the subspace J-'T is orthogonal (in 
the sense of the duality product) to the space of transmitted wrenches CW, 
that is 

FT = CW^ (25) 

We have defined the spaces of freedom twists and constraint wrenches as 
subspaces of of the Lie algebra se(3) and its dual. However these spaces 
express constraints on the twists and wrenches of the rigid bodies related 
by the kinematic pairs and hence are expressed in some common frame with 
configuration Qkp- (In most cases equal to the configuration of one of the 
related bodies. ) Consequently, the constitutive relations of a kinematic 
pair is given in terms of its pair of twists and wrenches {Txp, Wkp) G 
Tq^pSE (3) X T^^pSE (3) as follows 

Ad*Q^^ Wkp G CW and Ad^-i Tkp G FT (26) 

Hence the constitutive equations of a kinematic pair may be expressed as 
the following non-constant separable Dirac structure: 

-Dew {Qkp) = {{Tkp, Wkp) G Tq^^SE (3) x T*q^^^SE (3) , 

Ad*n Wkp G CW and Ad^-i Tkp G CW^ \ 



The kinematic pair introduced above represents ideal kinematic constraints; 
in general, however, mechanical work may be produced at the kinematic pair 
due to the presence of actuators or springs and dampers. Such an interaction 
is captured by considering the linear space IW := se*(3)/CW (which may be 
identified with a subspace of se*(3) complementary to the space of constraint 
wrenches CW). The space of interaction twists is then defined as its dual 
space XT := XW* ~ CW . Using the canonical projection vr of se*(3) onto 
IW, together with its adjoint map vr*, one may thus define an additional pair 
of port variables enabling to connect actuators, damper or spring elements 
to the kinematic pairs. 

The resulting interacting kinematic pair is then defined as a 2-port ele- 
ment with constitutive relations defined by the following non-constant sep- 
arable Dirac structure 

T^cw (Qkp) = {{Tkp, Wkp, Tj, Wi) G 
Tq,,SE{3) X T^^^SEiS) X CW^ x se*{3)/CW, (28) 

Wi = Tro Ad*Q^^ (Wkp) and Tkp = -Adg^.^ o tt* (T/)} 

It is easy to check that if Wj = then the interacting kinematic pair reduces 
to the kinematic pair as defined before. 

3.3.3 The kinestatic connection network 

The primary graph of the mechanism together with the kinematic pairs 
is called the kinestatic model of the mechanical system. Its associated Dirac 
structure is the composition of the Dirac structures corresponding to the 
kinematic pairs with the flow-continuouq^j graph Dirac structure of the pri- 
mary graph. 

Consider a mechanism defined by its primary graph Q composed of urb 
internal vertices (associated with the rigid bodies), Ub boundary vertices 
corresponding to rigid bodies with zero inertia operator and ukp edges (as- 
sociated with the kinematic pairs). Define the vertex space Aq 3 T , and 
the edge space Ai 3 T^^ with respect to the Lie algebra se (3), which rep- 
resent respectively the external twist of the rigid bodies and the kinematic 
pairs. The dual spaces A^ 9 W^^ , respectively A^ 3 W^^, then represent 
the external wrenches of the rigid bodies, respectively the wrenches of the 



^°Or the effort-continuous grapti Dirac structure in case ttie rigid bodies corresponding 
to ttie boundary vertices tiave non-zero inertia operator. 



kinematic pairs; see also Remark 2.1. The twists and wrenches of the bound- 
ary vertices (the rigid bodies with zero inertia operator) are associated with 
the vertex space A^ 3 T , respectively its dual A 9 W . Kirchhoff's laws 
on the twists and wrenches [lOj amount to constraining these variables to 
belong to the flow-continuous graph Dirac structure, i.e., 



(t^p, ly^^, T«^, T^^^, ^^ w") G VfiG) 



Composition of T>f{G) with the Dirac structures T>cy\;{QKp) corresponding 
to all the kinematic pairs then results in the Dirac structure "Dks of the 
kinestatic model: 

(t^, W^, T^^, W^^, T\ W'') G Vks (29) 

3.3.4 The dynamics of a mechanism 

The state space of the complete mechanism is the product space of the state 
spaces of all the rigid bodies, i.e., 

X = (5^(3) X se*(3))"«^ 3 x, (30) 

where ns denotes the number of rigid bodies (equal to the number of internal 
vertices of the primary graph) . 

Recalling that the rigid body dynamics is defined as a port Hamiltonian 
system with respect to the Dirac structure (|24p one then obtains the overall 
Dirac structure Pm of the mechanism by composing the Dirac structure 
"Dxs of the kinestatic model with the Dirac structures P_r_b of all the rigid 
bodies. Finally, defining the Hamiltonian Hm (x) as the sum of the Hamil- 
tonians of each body one obtains the following port-Hamiltonian model of 
the mechanism: 

-^,^{x),T\w\T\w'^eVM (31) 

3.4 Port-Hamiltonian formulation of consensus algorithms 

While all previous examples of port-Hamiltonian dynamics on graphs arise 
from physical modeling the system treated in this subsection has a different 
origin. Nevertheless, it shares the same structure, and in fact, turns out to 
have dynamics equal to the mass-damper system treated before. 

Consider a network of N agents moving in linear space 7^, whose interac- 
tion topology is described by an undirected graph Q (symmetric interaction). 



Denote by E{Q) the edges of this undirected graph, consisting of unordered 
pairs (f, It;) of vertices t',t(;. Hence {v,w) G E{Q) if and only if (w,t;) € E{Q). 
Thus the vertices of the graph correspond to the agents, and the edges to the 
symmetric interactions between them. Distinguish between leader and fol- 
lower agents, see e.g. [27], and associate the leader agents to the boundary 
vertices and the follower agents to the internal vertices. 

Associated to each agent v there is a vector x^ gTZ describing the motion 
in the linear space TZ. In the standard consensus algorithm, see e.g. [23] . 
the vector Xy of each follower agent v satisfies the following dynamics 

Xv{t) = - ^ s'(„^^)(x„(t) -x^(t)) (32) 

{v,w)£E{g) 

where fl'(^^^) > denotes a certain positive-definite weight matrix associated 
to each edge. 

For simplicity of exposition let us take the linear space TZ to be equal 
to M in the rest of this section, implying that 5f(t,^^) > are just positive 
constants. Collecting all follower variables Xy into one vector x G M % and 
all leader variables x^ into one vector u G M^'' , it is readily checked that the 
dynamics can be written as 

X = -BiGBjx - BiGBlu (33) 

with B the incidence matrix of the graph endowed with an arbitrary orienta- 
tion, and G the diagonal matrix with elements 5(t,^„,) corresponding to each 
edge (f , w). This defines a port-Hamiltonian system with respect to the flow- 
continuous graph Dirac structure T>f{Q). Indeed, consider the Hamiltonian 
function given by H{x) := ^ || x |p then p3]) is equal to 

X = -BiGBf — {x) - BiGBlu (34) 

which are the same equations as for the mass-damper system (j2ip . with 
u = Cb G A^ y = /fe G Afo, and ||(x) = x. 

Using the properties of the Laplacian matrix, see the discussion at the 
end of Section 13.2. H it can be shown [39l |33] that if the graph Q is connected 
then for any u there exists a unique x satisfying BiGBjx + BiGB^u = 
0, that is, an equilibrium of the system. Furthermore, this equilibrium is 
asymptotically stable. 

Note that the corresponding artiflcial output vector 

BB 

y := BbGBf — {x) + B^GBfu 



equals minus the rate of the leader variables if the leader variables were 
supposed to obey the consensus algorithm with regard to the follower agents 
(which is not the case). Hence this artificial output measures the discrepancy 
between the leaders and the followers. 

3.5 Port-Hamiltonian formulation of coordination control 

Another example of port-Hamiltonian systems on a graph is provided by the 
coordination control of A^ dynamical systems |^. Consider an open graph 
where each internal vertex corresponds to a port-Hamiltonian system. Coor- 
dination will be sought by designing a port-Hamiltonian dynamics associated 
to each edge of the graph. 

For clarity of exposition we will restrict attention to the flow- continuous 
graph Dirac structure; the treatment for the effort-continuous graph Dirac 
structure runs analogously. Thus we consider an open graph G, where to each 
internal vertex Vi E Vj we associate a Dirac structure P^. , a Hamiltonian 
Hy.{xy.), and a resistive relation Ry.. That is, for each internal vertex Vi we 
have a port-Hamiltonian system 

with vectors of external flows and efforts fy. , Cy. , where the flows and efforts 
f^, Cy. corresponding to the resistive elements are related by the dissipative 
relation 

RvACe^) = (36) 

All these port-Hamiltonian systems are coupled to the flow-continuous graph 
Dirac structure T>f{Q) by the standard interconnection 

fv, = -foi;v,, e„, = eyl, Vi G Vi (37) 

where foi-Vi and e|]* are the components of the vector /oj, respectively e^^, 
corresponding to the internal vertex Vi. What results is a port-Hamiltonian 
system, with Dirac structure being the composition of the internal vertex 
Dirac structures Dy- with the flow-continuous graph Dirac structure Vf^Q), 
still having external flow and effort variables 



^'^The set-up is heavily inspired by the innovative paper [I]. A major difference is that 
in our case the dynamics corresponding to the vertices and edges are all assumed to be 



port-Hamiltonian, instead of merely passive as in jT]. 



Coordination control is now sought by interconnecting the edge flow and 
effort variables /i,e^ to port-Hamiltonian systems attached to the edges. 
The simplest possibility for doing this is to associate energy storage to the 
edges, that is 

/i = -Xe, e^ = -^— ^(xe), e G f , (38) 

resulting in a port-Hamiltonian system with total Hamiltonian 

Htot = Yl H-^(^-^) + Yl ^^(^^) ^39) 

The aim is now to design H^, in such a manner that the minimum of Hi^t 
corresponds to a desired set-point. Then under the additional assumption 
of enough dissipation in the system the system will converge to this desired 
configuration; see [1] for further details in a passivity framework. 

4 Port-Hamiltonian systems on graphs obtained 
from symmetry reduction of symplectic Hamil- 
tonian systems 

In this section we will show how the port-Hamiltonian formulation of the 
mass-spring system on an open graph can be alternatively obtained from 
a canonical symplectic formulation by symmetry reduction, exploiting the 
invariance of the Hamiltonian function (in particular, of the spring potential 
energies) . 

4.1 Symmetry reduction from the symplectic formulation 



Let us return to the formulation of a mass-spring system in Section 
where the vertices correspond to the masses, and the edges to the springs 
in between them. An alternative is to consider the configuratioro vector 
Qc G A*^ =: Qc, describing the positions of all the masses. In fact, this is 
the classical starting point for Lagrangian mechanics, where we do not start 
with the energy variables q and p, but instead we start with the configuration 
vector Qc and the corresponding velocity vector Qc- The classical Hamiltonian 
formulation is then obtained by defining the vector of momenta p G Aq = Q* 



^^Note that qc G A'' is defined to be a function gc : V — >■ R, assigning to each vertex its 
configuration in 7?.. 



as p = Mqc (with M the diagonal mass matrix), resulting in the symplectic 
phase space 

QcxQ*=A^x Ao 

For ordinary springs the relation between q^ € A'' and the vector q £ A^ 
describing the elongations of the springs is given as q = B qc- Hence in this 
case the Hamiltonian can be also expressed as a function He of {qc,p) by 
defining 

H,{q„p):=H{B^q„p) (40) 

It follows that the equations of motion of the mass-spring system (with 
boundary masses) are given by the canonical Hamiltonian equations 



^fe,p) 




I 

s' = [0 l]^(Qc,i)' 



dHc 
dp 



h (41) 



where, as before, fh are the external forces exerted on the boundary masses 
and Cfe are their velocities. 

What is the relation with the port-Hamiltonian formulation given in 
Section [3?2]? It turns out that this relation is precisely given by the standard 
procedure of symmetry reduction of a Hamiltonian systeno- Indeed, since 
B \ = the Hamiltonian function Hc{qc,p) given in (I40p is invariant under 
the action of the group = R acting on the phase space A*^ x Aq = M^^ by 
the symplectic group action 

{qc,p)^{qc + al,p), a G 6 = M (42) 



From standard reduction theory, see e.g. [18j and the references quoted 
therein, it follows that we may factor out the configuration space Qc '■= A^ 
to the reduced configuration space 

Q := AV6 (43) 



^^This relation can be regarded as the discrete, graph-theoretic, version, of the corre- 
spondence between the port-Hamiltonian formulation of the Maxwell equations (using the 
Stokes-Dirac structure) and its symplectic formulation using the vector potential of the 
magnetic field, cf. [T51H0] . 



Let us first assume that the graph is connected, in which case we have the 
equality, cf . ()54p , ker B^ = span 1 . Then we have the following identification 

Q := A7© ~ 5^A° C A^ (44) 

Hence the reduced state space of the mass-spring system is given by im B^ x 
Aq, where im B C A^ . (Note that im B = A^ if and only if the graph does 
not contain cycles.) Furthermore, under the symmetry action the canonical 
Hamiltonian equations (|^T]) on the symplectic space A" x Aq reduce to the 
port-Hamiltonian equations (j45p on im B^ x Aq C A^ x Aq obtained before: 



B^q, = B^^{q,,p)=BTmq,p) 



dp wCiF/ Qp 



p = -^fe,p) + 

e' = [0 ^]f(?'P) 



fb = -B^iq,p) + 



dq 



h (45) 



dp 

In case the graph is not connected, then the above symmetry reduction can 
be performed for each component of the graph (i.e., the symmetry group is 
M^, with kg denoting the number of components of the graph Q), yielding 
again the reduced state space miB^ x Aq. 

Since the port-Hamiltonian equations (I45p for a connected graph can 
thus be regarded to be the result of a 1-dimensional symmetry reduction, it 
follows from general reduction theory that a conserved quantity is present. 
Indeed the total momentum ptot := l^P satisfies the balance law 

j^Ptot = llh (46) 

In particular, ptot is conserved if the sum of the external forces l^/h is equal 
to zero. 

Remark 4.1. In addition to the total momentum there may be other con- 
served quantities. Indeed, considering the Poisson structure matrix defined 
in liTp it follows that all quantities q with Bq = (corresponding to the 
cycles of the graph Q) are conserved. 

5 The KirchhofF-Dirac structure on graphs and its 
port-Hamiltonian dynamics 

In this section we consider a third canonical Dirac structure on a graph, 
which corresponds to constraining the flows at the internal vertices to zero. 
Thus for the resulting port-Hamiltonian systems there is no energy-storage 
or dissipation associated with the vertices. 



5.1 The Kirchhoff-Dirac structure 

The Kirchhoff-Dirac structure is defined as 

^KiQ) := {{fi,e\fb,e'') E Ai x A^ x A^ x A^ | 
B^h = 0, Bhfi = fb, 3e0* G A^^ s.t. e^ = -B*e°' - B^e''} 



(47) 



Note that, in contrast to the flow/effort-continuous graph Dirac structures, 
the Kirchhoff-Dirac structure only involves the flow and effort variables of 
the edge and boundary vertex spaces (not of the internal vertex spaces) . 

Proposition 5.1. T>k{Q) is a separable Dirac structure. 

Proof. The Kirchhoff-Dirac structure is equal to the composition of the 
fiow-continuoua^^l graph Dirac structure Df{Q) with the following trivial 
separable Dirac structure defined on the internal vertex spaces 

{{foi,e^')eAo^xA^'\foi = 0} 

The result then follows from Proposition 12.81 ■ 

Port-Hamiltonian dynamics with respect to the Kirchhoff-Dirac struc- 
ture is defined completely similar as in the case of the flow/effort-continuous 
graph Dirac structure; the difference being that energy-storing or dissipative 
relations are now only defined for the flow and effort variables corresponding 
to the edges. 

5.2 Electrical circuits 

The example of a port-Hamiltonian systerq^ with respect to a Kirchhoff- 
Dirac structure is an electrical RLC-circuit, with circuit graph Q. In this case 
the elements of Ai and A^ denote the vectors of currents through, respec- 
tively the voltages across, the edges. The Kirchhoff-Dirac structure in this 
case exactly amounts to the Kirchhoff's current and voltage laws (whence 
its name). Furthermore, the effort variables e^ are the potentials at the ver- 
tices, while the boundary flows and efforts /fe,e* are the boundary currents, 



^*0r the composition of the effort-continuous graph Dirac structure with {(/o,e'^) £ 
Ao X A° I /o = 0}. 

^^ Strictly speaking, the terminology 'port-Hamiltonian' is slightly confusing in this con- 
text, because 'ports' in electrical circuits are usually defined by pairs of terminals, that 
is pairs of boundary vertices with external variables being the currents through and the 
voltages across an edge corresponding to each such port. See also Remark [2.21 



respectively boundary potentials at the boundary vertices (the terminals of 
the electrical circuit). 

On top of Kirchhoff 's laws, the dynamics is defined by the energy-storage 
relations corresponding to either capacitors or inductors, and dissipative re- 
lations corresponding to resistors. The energy-storing relations for a capac- 
itor at edge e are given by 

Qe = -Ie, Ve = ^iQe) (48) 

with Qe the charge, and HceiQe) denoting the electric energy stored in 
the capacitor. Alternatively, in the case of an inductor one specifies the 
magnetic energy Hiei^e), where $e is the magnetic flux linkage, together 
with the dynamic relations 

ie = V;, -/, = ^($,) (49) 

Finally, a resistor at edge e corresponds to a static relation between the 
current /g through and the voltage Ve across this edge, such that Vele < 0. 
In particular, a linear (ohmic) resistor at edge e is specified by a relation 

Ve = -Rele, with Re > 0. 

Alternatively, we can decompose the circuit graph Q as the interconnec- 
tion of a graph corresponding to the capacitors, a graph corresponding to 
the inductors, and a graph corresponding to the resistors. For simplicity let 
us restrict ourselves to the case of an LC-circuit without boundary vertices. 
Define V as the set of all vertices that are adjacent to at least one capacitor 
as well as to at least one inductor. Then split the circuit graph into an open 
circuit graph Q corresponding to the capacitors and an open circuit graph 
gL corresponding to the inductors, both with set of boundary vertices V. 
Denote the incidence matrices of these two circuit graphs by 



B 



c 






B^ 



Br 

Bt 



Assuming for simplicity that all capacitors and inductors are linear we arrive 
at the following equations for the C-circuit 

B^Q = I^, BfQ = 0, 






BGT^C ^ C-^Q-Bf^ijf 



with Q the vector of charges of the capacitors and C the diagonal matrix 
with diagonal elements given by the capacitances of the capacitors. Similarly 



for the L-circuit we obtain the equations 



$ = 


- Bt^^pi: + B^i;^ 


= 


= Bl^L-^^ 


tL _ 

b - 


= -Bj^L-^^ 



with $ the vector of fluxes and L the diagonal matrix of inductances of all 
the inductors. 

The equations of the LC-circuit are obtained by imposing the inter- 



connection constraints ip^ = il)^ =: t/ji and If^ + Ig' = 0. By eliminating 



the boundary currents If^,I^ one thus arrives at the differential-algebraic 
port-Hamiltonian equations 



c 



B' 



B^ 





B^ 
Bt. 



-6 



-Q 




BCT- 






For a formulation of pure R,L or C circuits, and their weighted Laplacian 
matrices, we refer to |33] . 



5.3 Mass-spring systems with regard to a Lagrangian tree 

An alternative port-Hamiltonian formulation of mass-spring(-damper) sys- 
tems, in terms of the Kirchhoff-Dirac structure, can be given as follows. 
Recall the port-Hamiltonian formulation on A^ x Aq with respect to the 
effort-continuous graph Dirac structure ^^{Q), in which case the masses 
correspond to the vertices, and the springs to the edges of the graph Q, 
which we assume to be connecteqij. This graph can be extended to an 
augmented graph ^aug by adding a ground vertex g and adding edges from 
every vertex v ol Q towards this ground vertex. (The augmented graph is 
called a Lagrangian tree.) Furthermore, by constraining the effort eg at the 
ground vertex to be zero we can equate the efforts e„ at every vertex v of G 
with the effort eyg at the edge from u to g' of the augmented graph ^aug- In 
this way we can identify the effort-continuous graph Dirac structure T>e{Q) 
with the Kirchhoff-Dirac structure 'DK{Qimg) with the additional constraint 



For non-connected graphs Q the same construction can be done for every connected 
component. 



eg = 0. (Note that this is again a separable Dirac structure since it equals 
the composition of the Kirchhoff-Dirac structure T>K{Ga,ug) with the trivial 
Dirac structure {{fg,eg) \ eg = 0}.) 

In this way, the masses become associated with the edges eyg from every 
vertex v to the ground vertex g. The interpretation of the ground vertex g is 
that it represents the reference point (with velocity eg being zero). The flow 
fg at the ground vertex g equals the total force exerted on a mass located 
at this reference point. 

5.4 Hydraulic networks 

The interpretation of the flow-/effort-continuous graph Dirac structures and 
the Kirchhoff-Dirac structure as capturing the basic conservation/balance 
laws of a network becomes especially tangible for hydraulic networks. 

A hydraulic network can be modeled as a directed graph with edges 
corresponding to pipes, see e.g. [28l|TT]. The vertices may either correspond 
to connection points with fluid reservoirs (buffers), or merely to connection 
points of the pipes. Let Xy be the stored fluid at vertex v and let qe be the 
fluid flow through edge e. Collecting all stored fluids x^ into one vector x, 
and all fluid flows q^ into one vector q, the mass-balance is summarized as 

X = -Bq (50) 

with B denoting the incidence matrix of the graph. In the absence of fluid 
reservoirs this simply reduces to Kirchhoff 's current laws Bq = 0. 

For incompressible fluids a standard model of the fluid flow qe through 
pipe e is 

Jeqe = P^- Pj " KiQe) (51) 

where Pj and Pj are the pressures at the tail, respectively head, vertices of 
edge e. Note that this encapsulates in fact two effects; one corresponding 
to energy storage and one corresponding to energy dissipation. Defining 
the energy variable ipe := JeQe the stored energy in the pipe associated 
with edge e is given as 2T^e = ^JeQe- Secondly, Xe{qe) is a damping force 
corresponding to energy dissipation. 

In the case of fluid reservoirs at the vertices the pressures P^ at each ver- 
tex V are functions of Xy, and thus, being scalar functions, always derivable 
from an energy function 

P, = ^{x,), f G V, (52) 

OX,, 



for some Hamiltonian Hy{xy) (e.g. gravitational energy). The resulting dy- 
namics (with state variables x and q) is port-Hamiltonian with respect to the 
graph Dirac structure T>f{Q) = T>e{Q). The set-up is immediately extended 
to boundary vertices (either corresponding to controlled fluid reservoirs or 
direct in/ or out-flows). 

In the absence of fluid reservoirs {Bq = 0) the dynamics reduces to 
a port-Hamiltonian system in the state variables q with respect to the 
Kirchhoff-Dirac structure . 

Remark 5.2. Similar models are used in inventory theory and supply-chain 
applications. 

5.5 Properties of the boundary flows and efl'orts of the Kirchhoff- 
Dirac structure 

The fact that the internal vertex flows in the definition of the Kirchhoff-Dirac 
structure are all zero (and consequently no storage or dissipation at the 
vertices takes place) has a number of specific consequences for the behavior 
of the boundary flows and efforts (see [12] for closely related considerations). 
Assume (for simplicity of exposition) that 7^ = M. It is a well-known 
property [5] of any incidence matrix B that 

I'^B = (53) 

where 1 denotes the vector with all components equal to 1. Furthermore [5], 
the rank of the incidence matrix is equal to the number of vertices minus 
the number of components P^l of the graph Q. In fact, each component of 
the graph satisfies the property 

ker S^ = span 1 (54) 

with B the (restricted) incidence matrix of this component and the dimen- 
sion of 1 equal to its number of vertices. This implies several properties for 
the previously defined Dirac structures on the graph; we will concentrate on 
its consequences for the Kirchhoff-Dirac structure. From the definition of 
the Kirchhoff-Dirac structure it follows that 

= l^Bfi = llBbfi = -llh (55) 

with 1ft denoting the vector with all ones of dimension equal to the number of 
boundary vertices. Hence the boundary part of the Kirchhoff-Dirac structure 



^'^A component is a maximal subgraph which is connected, that is, every two vertices 
are hnked by a path of, -non-oriented-, edges. 



of an open graph is constrained by the fact that the boundary flows add up 
to zero. Dually, we may always add to the vector of vertex efforts e^ the 
vector 1 leaving invariant the edge efforts e^ = B e^ . Hence, to the vector 
of boundary efforts e^ we may always add the vector 1^. 

Proposition 5.3. Consider an open graph Q with Kirchhoff-Dirac structure 
Vk{G)- Then for each (/i,e^, /t, e*) G Vk{Q) it holds that 

tlfb = 0, 

while for any constant c € M 

{fi,e\fb,e' + clb)eVK{g) 

This proposition implies that we may restrict the dimension of the space 
of boundary flows and efforts A;, x A** of a connected graph by two. Indeed, 
we may define 

Afered := {fb G Afo I /ft G ker 1^} 

and its dual space 

A^,d:=(Afered)* = AVimlft 

It is straightforward to show that the Kirchhoff-Dirac structure VxiG) re- 
duces to a linear subspace of the reduced space Ai x A^ x Abred x ^'^ed' 
which is also a Dirac structure. An interpretation of this reduction is that 
we may consider one of the boundary vertices, say the first one, to be the 
reference vertex, and that we may reduce the vector of boundary efforts 
e* = (e''^ • • • , e**) to a vector of voltages (e*^ - e*\ • • • ,e^^-e^^). A graph- 
theoretical interpretation is that instead of the incidence matrix B we con- 
sider the restricted incidence matrix [3j. 

For a graph G with more than one connected component the above holds 
for each connected component. It follows that there are as many indepen- 
dent constraints on the boundary flows /{, as the number of the connected 
components of the open graph Q. Dually, the space of allowed boundary 
efforts e is invariant under translation by as many independent vectors If, 
as the number of connected components. 

A complementary view on Proposition 15.31 is the fact that we may close 
an open graph Q to a closed graph G as follows. Consider first the case that 
Q is connected. Then we may add one virtual ground vertex Vg, and virtual 
edges from this virtual vertex to every boundary vertex Vb € Ve, in such 



a manner that the Kirchhoff-Dirac structure of this graph Q extends the 
Kirchhoff-Dirac structure of the open graph Q. In fact, to the virtual vertex 
Vg we may associate an arbitrary potential e^{vg) (the ground-potential), 
and we may rewrite the externally supplied power < e^ | /;, > as (since by 

m E., h{vb) = 0) 

<e'\h >= Y.^e\v,) - e\vg)h{v,) = ^ e^\v^)Mv^), (56) 

where e^^{vi)) := e^{v!,) — e^{vg) and fib{vb) '■= fh{vh) denotes the effort 
across and the flow through the virtual edge towards the boundary vertex 
Vfe. It is clear that for every element (/i,e^, /fe, e'') € VxiG) corresponding 
to the open graph Q there exists e^^ such that (/i, e^, /ifc,e^^) G Q for the 
closed graph Q, and conversely for every {fi,e^, fib,e^^) G Q there exists 
e* such that (/i, e^, /;,, e*) G Vk{G)- This construction is extended to non- 
connected graphs by adding a ground vertex to each component containing 
boundary vertices. 

6 Port-Hamiltonian dynamics on complexes 

In this section the formulation of port-Hamiltonian dynamics on graphs will 
be extended to k- complexes. A main motivation for doing so is to be able 
to represent (continuous) systems of conservation laws directly on a discrete 
spatial domain by using smplicial comp/execj, in such a way as to retain 
the geometric nature of the involved variables. 

This is closely related to an intense research activity on the use of mixed 
finite elements and discrete differential forms for the spatial discretization of 
partial differential equations, see among many other references e.g. [6l I14j . 
The novelty of the approach taken in this section resides in systematically 
taking into account the boundary interaction of such systems and in the re- 
sulting formulation of the discrete conservation laws in terms of Dirac struc- 
tures and finite-dimensional port-Hamiltonian systems. More specifically, 
our approach can be regarded as the discrete analog of the formulation of 
systems of continuous conservation laws with open boundary as distributed- 
parameter (infinite-dimensional) port-Hamiltonian systems with respect to 
a special type of Dirac structure, the Stokes-Dirac structure, as performed 
in [36, il7J. In the absence of a boundary (or for fixed boundary conditions) 



^^Also other complexes, e.g. corresponding to cubes in R"^ instead of tetrahedra, could 
be considered. 



this is in its turn closely related to the extensive literature on the Hamilto- 
nian formulation of systems of evolution equations using infinite-dimensional 
Poisson brackets; see e.g. [H] and the references quoted therein. 

Apart from the main motivation of structure-preserving spatial discretiza- 
tion we will also discuss at the end of this section the extension of consensus 
algorithms to A;-complexes. 

6.1 Prom graphs to /c-complexes 

A directed graph with incidence matrix B can be regarded as an example 
of what is called a 1- complex. Indeed, the sequence 

Ai 4 Ao ^ M 

satisfies the property 1 o B = 0. More generally, a A;-complex A can be 
specified by a sequence of real linear spaces Aq, Ai, • • • , A^, together with a 
sequence of linear incidence operators 

Afc 4 Afc_i -^ • • • Ai 4 Ao 

with the property that dj-iodj = 0, j = 2,- ■ ■ ,k. The vector spaces Aj, j = 
0, 1 • • • ,k, are called the spaces of j- chains. Each Aj is generated by a finite 
set of j-cells (like edges and vertices for graphs) in the sense that Aj is the 
set of functions from the j-cells to M. 

A typical example of a A;-complex is the triangularization of a A;-dimensional 
manifold, with the j-cells, j = 0, 1, • • • ,k, being the sets of vertices, edges, 
faces, etc. (the so-called simplicial complex). 

Example 6.1. Consider the triangularization of a 2- dimensional sphere by 
a tetrahedron with 4 faces, 6 edges, and 4 vertices. The matrix representation 
of the incidence operator 82 (from the faces of the tetrahedron to its edges) 
is 

< ?;i?;2 > 1 -1 

< fifs > -1 1 

< ■Uif4 > -1 1 

< V2V3 > 1 0-1 

< V2V4 > -1 1 

< f3f4 > 1 0-1 



V4 



V3 




Vi 



Figure 1: Tetrahedron triangularizing a sphere 

where the expressions < wiW2f3 >,••• denote the faces (with corresponding 
orientation), and < wit;2 >, • • • are the edges. The matrix representation of 
the incidence operator di (from edges to vertices) is given as 

< W1U2 > < ViV^ > < V1V4 > < V2VS > < V2V4 > < V^Vi > 



<Vl> 


-1 


-1 


-1 











<V2> 


1 








-1 


-1 





<V3> 





1 





1 





-: 


< W4 > 








1 





1 


1 



It can be verified that 9i o O2 = 0. 

Denoting the dual Unear spaces by A^ , j = 0,1 
following dual sequence 



,k, we obtain the 



where the adjoint maps d^ j = 0, 1 • • • ,k, satisfy the analogous property 
d^ o d^^^ = 0, j = 2,- ■ ■ ,k. The elements of A-' are called j-cochains. 



6.2 Canonical Dirac structures associated with open k- 
complexes and port-Hamiltonian systems 

An open /c-complex is identified by identifying a subset V?^_^-) of the set 
of all {k — l)-cells, called the boundary {k — l)-cella^^l. while the remaining 
(k — l)-cells constitute the set V^ _-^n of internal {k — l)-cells. Define the 
linear space of functions from this subset of {k — l)-cells to M as A^ C Afc_i 
with dual space denoted as A''. Decompose correspondingly 9^ : A^ — > Afe_i 
as dk = {dk-i,dk-i,), with adjoint mapping d'' = {d''''^ , d'''^) . The space of 
flows and efforts {fb, e ) G A;, x A describes the distributed terminals of the 
open /c-complex. 

Similar to the definition of the flow-continuous graph Dirac structure, 
see (|10p . we define the flow- continuous k-complex Dirac structure as the 
subspace 



X 



P/(A) := {(/fc,e\/fc_i;i,e'=-i^/fe,e'') G Afc X A'^ X Afc_i;, X A'= 
A, X A'' I dk.ifk = fk-i;u dk-bfk = fb, e^ = -d'^^^e^-i'^ - d^'^e'} 

Similarly we can define the eff^ort- continuous k-complex Dirac structure 
Pe(A) as the subspace 

Pe(A) := {(/fc,e^/fc_l;i,e'=-l'^/fe,e*) € A^ x A'= x A,,_i x K^-^ x A^ x A'' | 

difk = fk-i, dk-bfk = fk-i;b + fb, e^ = -d^-^'e^-\e'' = e'^-^'''} 

(58) 

As in the graph case, Proposition l2.7l implies that Pj ( A) , Pg ( A) are separable 
Dirac structures. Furthermore, the linear subspace 

Pi^(A):={(/fc,e^/^„e^)GAfcxA'^xAfexA^| 

(59) 
^k■,^fk = 0, dk;bfk = fb, ae'^-i'^ s.t. e^ = -d'^'^e^-^'^ + d'^'^e^} 

defines a separable Dirac structure, called the Kirchhoff- Dirac structure on 
k-complexes] the proof is similar to Proposition 15.11 Analogously to the 
graph case, the Kirchhoff-Dirac structure for open A;-complexes implies cer- 
tain constraints on the incoming 'currents' fb- Indeed, by the fact that 
9fc_i o dk = it follows that d(^k_iy^bfb = 0, where 5(fc_i).f, denotes the 
{k — l)-th incidence operator restricted to A^ C Afc_i. (Note that in the 



^^One could also consider as boundary cells subsets of the j'-th cells for j ^ fc — 1. 
In particular, choosing j — k would correspond to 'distributed interaction'. The choice 
j — k ~ 1 corresponds to the important case of 'boundary' interaction. 



case of a graph the role of d(^k_iyi, is played by the linear map 1^.) As in 
the case of graphs, this allows us to reduce the Kirchhoff behavior to a space 
that is still a Dirac structure, or, alternatively, to close the open A;-complex. 
This is done by completing the open /c-complex A with space of boundary 
currents A^ by an additional set of {k — l)-cells and /c-cells. 

Similarly to open graphs , we can define interconnections of open com- 
plexes. The most common form of interconnection of two open fc-complexes 
is through shared boundary (k — l)-cells, similarly to the interconnection of 
open graphs through shared boundary vertices o. 'Shared' here means that 
the {k — l)-subcomplexes of both open /c-complexes defined by the boundary 
(k — l)-cells are the same. 

Similarly to the graph case, port-Hamiltonian dynamics on a fc-complex 
can be defined in various ways; either with respect to its flow/effort-continuous 
/c-complex Dirac structure, or with respect to its Kirchhoff-Dirac structure. 
The next sections will describe a couple of examples. 

6.3 2D transverse magnetic Maxwell's equations defined on 
a 2-coniplex 

Consider the transverse magnetic Maxwell equations with respect to some 
plane in physical space. Such a situation arises when modeling a thin 
isotropic dielectric film comprised between two planar conductors, while 
neglecting boundary effects at the boundaries of the conducting plates. The 
electrical field may then be approximated as being orthogonal to the plane 
and the magnetic field as being tangential to the plane. Various discretiza- 
tion methods have been used for these transverse Maxwell's equations. Here 
instead we shall directly obtain a discretized model expressed as a port- 
Hamiltonian system defined with respect to the flow-continuous 2-complex 
Dirac structure. 

First we start with the simplicial 2-complex defined by a triangulariza- 
tion of the plane. Using the planar symmetry of the field, we shall associate 
different chains and cochains with the electrical and magnetic fields. 

Because of the assumption of orthogonality of the electrial field with the 
plane, we identify the electric field induction with a 2-cochain D G A^ and 
the electric field intensity with a 2- chain i? € A2. Furthermore, using the 
assumption that the magnetic field is tangential to the plane, we identify 



^°In principle there are several possibilities for the interconnection of fc-complexes, de- 
pending also on the definition of the boundary cells. For example, one may also define the 
interconnection of a fc-complex with an /-complex through a number of shared boundary 
p-cells, with p < k,p < I. 



the magnetic field induction with a 1- chain i? € Ai and the magnetic field 
intensity with a 1-cochain H & A^. 

Maxwell's equations applied to the discrete fields defined above may then 
be expressed, in the absence of a boundary, as the following two relations 
(conservation/balance laws) 

^ = d^H 

(60) 
^ = -d2E 

The properties of the dielectric is defined by the positive-definite electric 
permittivity tensor e on the space of electrical intensity 2-chains and by 
the positive-definite magnetic permeability tensor /x on the space of the 
magnetic field intensity 1-cochains. With an abuse of notation the map 
induced by the permittivity tensor from the electric intensity 2-chains to 
the electric induction 2-cochains will be also denoted by e. In the same way 
we use /i also for the induced map from the magnetic intensity 1-cochains 
to the magnetic induction 1-chains. This leads to the following constitutive 
relations 

D = eE B = fiH (61) 

and inversely 

E = e'^D H = fi-^B (62) 

Note that the discrete electric permittivity and magnetic permeability ten- 
sors depend not only on the continuous properties of the material but also on 
the geometry of the meshing (triangularization) . Furthermore, these tensors 
define the electromagnetic energy as (we use IH for the Hamiltonian function 
in order to distinguish it from the magnetic field intensity H) 

MiD, B) = ^ {e-'D, Z?> + i {,i-'B, B) (63) 

where ( , ) denotes the duality pairing between dual vector spaces. Note 
that E = ^{D,B) and H = ^{D,B). 

Introducing boundary edges (1-cells) we are then led to the following 
discrete Maxwell's equations 

dD ^ d^■,^m^D^B) + d''■''>e^ 

^ = -d2;^{D,B) (64) 

fb = d2;b^{D,B) 



which defines a port-Hamiltonian system with respect to the flow-continuous 
2-complex Dirac structurq^J, and with boundary port variables fh and e^ 
being respectively the magnetic flux and magnetic field intensity at the 
boundary edges of the system. 

6.4 Diffusive systems on A;-complexes 

In this subsection we will consider a particular type of port-Hamiltonian 
dynamics based on the fiow-continuous A;-complex Dirac structure, by spec- 
ifying dissipative relations between Afc_i.j and A'^"^'*, together with one 
type of energy-storage relations between A^ and A . The resulting class of 
systems will be called diffusive systems. 

Consider a A;-complex, together with its fiow-continuous Dirac structure 
Pa. We associate to every k-cell an energy storage, by imposing the following 
dynamical energy-storing relation between /^ and e^ 



„k 



dH 



—X = e 



f, = —{x), xEA'^ (65) 



dx 



-d^^''R^k,^^{x)-d^■^^e^ 



with H{x) the total stored energy, and a: G A the total vector of energy 
variables. Furthermore, to every internal {k — l)-cell we associate a (linear) 
resistive relation, leading to 

gfc-i;i ^ -Rf^_^..^ R = R^>0 (66) 

Then we arrive at the dynamics 

dH 

This defines a port-Hamiltonian system with inputs e'' and outputs /fe, which 
can be regarded as a standard model for a discretized diffusive system. For 
the continuous analog of this diffusive system, using differential forms in- 
stead of (co-)chains, we refer to [30] , 

As a typical example consider the discrete formulation of a mass-diffusion 
process on a bounded 3-dimensional spatial domain Z with boundary dZ 
(e.g., the diffusion of a chemical species in a medium due to the gradient of its 
chemical potential [U chap. 17]); see [37] for the example of discretized heat 



^^Note that, in order to stick to the same sign convention as in the continuous Maxwell's 
equations, the sign convention for the flow-continuous 2-complex Dirac structure is just 
the opposite to the one used above; see ([57|) . 



flow. Discretize the spatial domain Z by triangularization, corresponding to 
a standard simplicial complex consisting of the sequence of cells : 

A3 ^ A2 ^ A, ^ Ao ^M (68) 

tctrahcdra faces edges vertices 

To every tetrahedron we may now associate a variable denoting the total 
mass (of the chemical species) in this tetrahedron, leading to a vector x G A^ . 
Bookkeeping of the ingoing and outgoing mass flow in every tetrahedron 
corresponds to the mass balance equation (a typical example of a discrete 
conservation law) 

X = d^-^'N' + d^'^''N\ 

where iV* G A^* denotes the mass flux of the chemical species through the 
internal faces and N G A^ is the mass flux through the boundary faces 
(corresponding to a triangularization of the 2-dimensional boundary dZ). 
Associated to x there is the Hamiltonian H{x) (the free Gibbs energy), 
defining the vector /i = ^(x) G A3 of chemical potentials, which is conju- 
gate to the vector x G A^. Furthermore, the expressions 

Fi = <93;i/U G A2;i, Fh = ds.,b^J' G A2;b 

define the generating functions (resulting from non-equilibrium) at the inter- 
nal, respectively boundary, faces. Finally the internal generating functions 
determine the internal mass flux A^*, by a dissipative relation of the form 

A^^ = -RFi 

for a certain matrix R = R^ > 0. All this leads to the finite-dimensional 
port-Hamiltonian system 

X = -d^'iRd^.i^{x) + d^'^N'' 

^u (69) 

Fb = 53;b||(x) 

modeling a diffusive system with input N^ (mass flux through the boundary 
faces), and output Ff, (generating function at the boundary faces). 

6.5 Generalized consensus algorithms on fc-complexes 

Note that the above discretized diffusion model is in some sense opposite 
to the formulation of the standard consensus algorithm on a graph (or, 



equivalently, the model of a mass-damper system) treated before. Indeed, 
in the latter case, energy-storage is associated with the vertices (0-cells), 
while dissipation is associated with the edges (1-cells). This indicates that 
there are different possibilities to extend consensus algorithms from graphs 
to a higher-dimensional situation. Indeed, we may as well consider the 
generalized consensus algorithm corresponding to a diffusive system without 
boundary 



-dkRdkX, X G a'' ~ Aa 



(70) 



with R a positive diagonal matrix of appropriate dimensions, or, alterna- 
tively. 



z = -dkGdkZ, z € Afe_i ~ A' 



fc-i 



(71) 



with G a positive diagonal matrix. The dynamical properties of both pos- 
sibilities are rather different. For example, in the case k = 2 the state 
X G A^ ~ A2 for the first algorithm will converge to the subspace ker92) 
while for the second algorithm z € Ai ~ A^ will converge to the subspace 
ker d2 which contains the whole subspace imdi (and will be equal to it if 
the Betti numbers are zero). 

As an example, consider the triangularization of a 2-sphere by a single 
tetrahedron, with incidence structure given before. Compute d2Gd2 as 
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the generalized consensus algorithm i 
subspace ker ^2 = im di . 



(72) 
-d2Gd2Z converges to the 3-dimensional 



On the other hand, the generalized consensus algorithm corresponding 
to the diffusive system x = —d2Rd2X is given aa^i 
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(73) 



which converges to the 1-dimensional subspace spanl in M^. 

7 Conclusions 

We have laid down a general geometric framework for the description of 
physical networks dynamics on graphs. Starting point are the conservation 
laws corresponding to the incidence matrix of the graph. This defines three 
canonical Dirac structures on the combined vertex, edge, and boundary 
spaces and their duals, where the last one (the Kirchhoff-Dirac structure) 
corresponds to the absence of energy storage or energy dissipation at the 
vertices. Relating the internal flows and efforts by either energy-storing 
or dissipative relations yields various forms of port-Hamiltonian dynamics. 
We have illustrated the approach on a number of typical physical examples. 
Other examples that have not been discussed include e.g. supply-chain mod- 
els and compartmental systems. Furthermore we have shown how examples 
from a different origin, such as consensus algorithms and coordination con- 
trol strategies, can be formulated and analyzed within the same framework. 
In a future paper we will explore the potential of the framework for the 
modeling and analysis of chemical reaction networks, which cannot directly 
be modeled by standard directed graphs. 

In the last part of the paper we have extended the framework from di- 
rected graphs to general /c-complexes. We have shown how this allows to 
a give a spatially discretized model of the 2-D Maxwell equations and of 
general diffusive systems. Pertinent questions concern the relation of these 



^^This dynamics turns out to be equal to the standard consensus algorithm on the 
complete graph with 4 vertices. This can be also seen by associating to each of the 
faces of the 2-complex a vertex, and by assigning an edge between two vertices if the 
corresponding faces are adjacent. Since in the case of a single tetrahedron all faces are 
adjacent, a complete graph results. This yields a correspondence between the operator 82 
and the transposed incidence matrix B^ of the complete graph. 



methods to structure-preserving spatial discretization methods for their de- 
scription by partial differential equations models. Furthermore, we have 
indicated how the discrete diffusive systems models may motivate alterna- 
tive consensus algorithms. 

For clarity of exposition we have only considered the basic building blocks 
of port-Hamiltonian systems and graphs. Indeed, because the interconnec- 
tion of port-Hamiltonian systems again defines a port-Hamiltonian system, 
the framework also covers heterogeneous and multi-scale situations, where 
several of the constructs considered in the present paper are connected to 
each other. 

The models treated in this paper all correspond to conservation/balance 
laws within a particular physical domain. Furthermore, the energy-balance 
of the system components can be seen to result from the underlying con- 
servation laws and the assumption of integrable constitutive relations for 
energy-storage. On the other hand, port-based (bond-graph) modeling as 
originating in the work of Paynter [26j is aimed at providing a unifying mod- 
eling framework for multi-physics systems, by directly starting from energy- 
flows between system components from different physical domains. This also 
results in port-Hamiltonian models as has been amply demonstrated in e.g. 
pm EU [35l [32l [l3] • It is well-known that bond-graph modeling involves an 
additional abstraction step (e.g., different electrical circuits may lead to the 
same bond-graph, and, conversely, different bond-graphs may correspond to 
the same electrical circuit). Furthermore, in the case of electrical circuits 
port-based modeling starts with a port description (pairs of terminals), in- 
stead of the more basic starting point of terminals corresponding to con- 
servation laws. Although in most situations the resulting port-Hamiltonian 
systems are the same this leaves some questions to be answered; see also 



Finally, a main motivation for physical system network modeling, apart 
from its importance in analysis and simulation, is its use for control and 
design purposes. Port-Hamiltonian systems theory has been successful in 
exploiting the physical structure for control and design purposes, see e.g. 
[321 I25j . using various forms of passivity-based control, control by intercon- 
nection, and tools originating in network synthesis theory. The applications 
of this control methodology to the broad range of port-Hamiltonian systems 
on graphs and /c-complexes as formulated in this paper is an interesting area 
for further research. 
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